Abstract. Star counts at high and intermediate galactic 
latitudes, in the visible and the near infrared, are used to 
determine the density law and the initial mass function 
of the thick disc population. The combination of shallow 
fields dominated by stars at the turnoff with deep fields 
allows the determination of the thick disc mass function 
in the mass range 0.2-0.8Mq. Star counts are compared 
with simulations of a synthesis population model. The fit 
is based on a maximum likelihood criterion. The best fit 
y—i ' model gives a scale height of 800 pc, a scale length of 
O , 2500 pc and a local density of 10 -3 stars pc -3 or 7.1 10~ 4 
^ ' M©pc~ 3 for Mv < 8. The IMF is found to follow a power 
law dN/dm oc m~ 05 . This is the first determination of the 
thick disc mass function. 
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1. Introduction 

Solving the problem of the origin of the thick disc of 
the Milky Way depends on accurate determination of its 
present characteristics. Overall analysis of density, kine- 
matics and chemical properties help the understanding of 
the physical processes involved. A first step was obtained 
with the finding that the thick disc was probably formed 
by a merging event on the thin disc early in the early age 



of the Milky Way ( riommcr-Larsen fc Antonuccio-Dclogu 



1993 ; Robin ct al., 1996| ). This hypothesis was motivated 
by the kinematic findings (no gradient, no discontinuity 
between thin disc and thick disc in rotation and velocity 
dispersion) and abundances, in particular the [O/Fe] and 
[Mg/Fe] ratios. These ratios implies a sudden decrease in 
star formation rate between the thick disc and thin disc 
formation, lasting for at least 1 Gyr but not more than 
3 Gyr ( Gratton et al., 2000 ). More work is required to 
measure an accurate density law and to obtain a detailed 
description of the stellar population, such as the initial 
mass function and age. 

The thick disc density law can reasonably be mod- 
eled by a double exponential, or a density law close to 
a sech 2 . Star counts are presently unable to distinguish 
between these hypotheses. Even the determination of the 
scale height and the local density causes difficulties due 
to a slight degeneracy between these two parameters, as 
shown in various published results. Different analyses have 
resulted in eit her high scale height and small local density 
(for example, Reid fc Majewski (1993 ): 1400 pc and local 
density of 2% of the disc) or smal l scale height and higher 
local density (Robin ct al. (1996): 760 pc and density of 
5.6%, [Buser et al. (1999) ): 910~pc and 5.9%), with several 
intermediate results. The number of observed stars is de- 
rived by the integral of the density law over the line of 
sight. For an exponential density law, the mean distance 
of the stars in a complete sample is roughly twice the scale 
height. The thick disc dominates star counts at distances 



between 2 and 5 kpc over the galactic plane. However, 
photometric counts are not accurate enough to estimate 
the distances of stars at the turnoff with an accuracy of 
even a factor of two. Moreover no accurate determination 
of the local density has ever directly been done, even with 
Hipparcos, because of the small proportion of the thick 
disc locally with regard to the thin disc. 

The determination of the initial mass function (IMF) 
of the thick disc is an important issue in the contro- 



versy about the universality of the IMF. 3calo (1998) and 
Kroupa (2001) find that it flattens at masses below 0.5M Q , 



Send offprint requests to: C. Reyle 

* Partly based on observations made at CFHT, and at the 
European Southern Observatory 



either in our Galaxy or in the LMC. Whereas Scalo points 
out the difficulty measuring an IMF and argues that the 
uncertainties on the determination could be of the order 
of the apparent variations of the IMF, Kroupa analyzes in 
detail the variations of the IMF slope and concludes that 
star formation in higher metallicity environments appears 
to produce relatively more low-mass stars, i.e. a steeper 
slope at low masses. If confirmed, one should expect the 
thick disc to have a shallower IMF slope than the thin disc 
on average. 

Until now no direct measurement of the thick disc IMF 
has been done. All along the main sequence, the thick disc 
population is easily distinguishable from the disc and the 
halo using a good temperature indicator like the V-I index, 
because at a given magnitude below the turnoff of the halo 
(~ V > 17), the blue side is dominated by the halo and the 
red side by the disc, with the thick disc in between. At high 
latitudes, thick disc stars become a sizeable population at 
magnitude about 14-15 in V in wide field star counts when 
a significant proportion of turnoff stars are detectable on 
the blue side of the colour distribution. One can reach 
the peak of the luminosity function (at about My=ll) at 
magnitude 21 for stars at 1 kpc. By combining shallow star 
counts dominated by turnoff stars with deep photometry 
one should be able to compute the IMF slope on a large 
mass range, between the turnoff at about m = 0.8 M Q and 
the maximum near m = 0.2 M Q . 

In this paper we address the problem of the thick disc 
density law together with its IMF. The two problems can- 
not be treated separately from star count analysis. We 
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use a large set of stellar samples (described in section | 



in the visible and the near infrared, at shallow and deep 



magnitudes, to investigate the thick disc luminosity func- 
tion, the local density and its scale height and scale length. 
This analysis has been feasible using a coherent model of 
population synthesis which takes into account the various 
photometric systems of the data, different basic hypothe- 
ses on the parameters to test, and allows us to disentangle 
the bias effects in star count samples (section ||). A max- 
imum likelihood test is used to estimate the thick disc 
parameters (section El). Results are given in section]^ and 
discussed in section g. 

2. Selected data sets 

Data sets at medium and high galactic latitudes have been 
selected. They combine shallow and deep star counts, in 
the visible and the near infrared. The main characteris- 
tics of visible data are summarized in table 1. The dots 
in figure [l] give their distribution in galactic coordinates. 
The shallow star counts are three wide fields from Ojha 
et al. (1996| , |1999| ), one field from phiu (198C| ) close to 
the galactic anticentre and one patch from the CNOC2 
Field Galaxy Redshift S urvey (|Yee et al., 200C|) The deep 
fields are one field from Borra fc Lepage (1986|), one field 



Bicnayme et al. (1987a[ ,|b|); |Haywood et al. (1997| ); [Robin 
et al. (200q )~ 



from Reid & Majewski (1993) close to the North Galactic 
Pole, six fie lds from the DMS ( [Hall et al., 1996| ; joT" 



et al., 1998) used to search for quasars, two fields from 
Bouvicr et al. (1998 ) dedicated to the study of Pleiades 
and Praesepe clusters, and another field near the North 
Galactic Pole (SA57) obtained at the CFHT by Creze et 
al. (in preparation). This field is the deepest one. It is 
complete and free from significant galaxy contamination 
up to V=24. 

We also used near infrared data in I and J bands from 
the D eep Near Infrared Sur v ey of the Southern Sky: DE- 
NIS dEpchtcin ct al., 1997| , |l999|) reduced at the Paris 
Data Analysis Center. DENIS fields, named strips, are 
12' in right ascension and 30° in declination. We have se- 
lected 26 strips distributed on the sky at latitudes greater 
than 30° , as shown by the lines in figure [|. At lower lat- 
itudes, the disc population becomes dominant compared 
to the thick disc population. Selected samples are portions 
of strips over which the density gradient in declination is 
negligible. They cover 1 to 3 square degrees and are com- 
plete up to 1=17.5 or 1=18, depending on the strip. 

The absolute visual magnitude of thick disc stars in 
the selected samples covers a wide range, from My ~ 4 
up to 13. It peaks at My ~ 4-5 in the DENIS fields, My ~ 
6-7 in the shallow fields, My ~ 9-10 in the deep fields and 
My ~ 12 in SA57. 

3. The model of population synthesis 

We have used a revised version of the Besangon model of 
population synthesis. Previous versions were described in 



The model is based on a semi-empirical approach, 
where physical constraints and current knowledge of the 
formation and evolution scenario of the Galaxy are used 
as a first approximation for the population synthesis. The 
model involves 4 populations (disc, thick disc, halo and 
bulge) each deserving a specific treatment. The bulge pop- 
ulation, which is irrelevant for this analysis, will be de- 
scribed elsewhere. 

3.1. The disc population 

A standard evolution model is used to produce the disc 
population, based on a set of usual parameters : an initial 
mass function (IMF), a star formation rate (SFR) and 



a set of evolutionary tracks (see Haywood et al., 1997, 
and references therein). The disc population is assumed 
to evolve over 10 Gyr. A set of IMF slopes and SFRs are 
tentatively assumed and tested against star counts. 

A revised IMF has been used in the present analysis, 
tuned with the most recent Hipparcos results: the age- 
velocity dispersion relation is from pomez et al. (19"97 ), 
the local luminosity function is from [lahreiss fc Wiclcn 
(1997) and an IMF power law dN/dm oc m~ a is adjusted 



to it, giving a slope a = 1.6 in the low mass range 0.08- 
0.5M©, in good agreement with [Mcra et al. (1996| ): a = 
2±0.5 and [Kroupa (2000a| ): a = 1.3±0.5. The scale height 
has been self-consistently computed using the potential 
obtained from the constraints on the local dynamical mass 
from Creze et al. (1998). This aspect will be described 



elsewhere (Reyle et al., in preparation). 

The evolutionary model fixes the distribution of stars 
in the space of intrinsic parameters : effective temperature, 
gravity, absolute magnitude, mass and age. These parame- 
ters are converted into colours in various systems through 
stellar atmosphere models corrected to fit empirical data 
(Lejeune et al., 1997, 1998). While some errors still remain 



in the resulting colours for some spectral types, the overall 
agreement is good in the major part of the HR diagram. 
For l ow mass stars in the near infrared, synthetic colours 
from Baraffe et al. (199£ ) have been used. 



Since the Haywood et al. (1997) model is based on evo- 
lutionary tracks at solar metallicities, inverse blanketing 
corrections are introduced to give to the disc a metal- 



licity distribution in agreement with the Twarog (1980) 



age/metallicity distribution (mean and dispersion about 
the mean). 



3.2. The thick disc population 

In the population synthesis process, the thick disc popula- 
tion is modeled as originating from a single epoch of star 



formation. We use Bergbush & VandenBerg (1992) oxy- 



gen enhanced evolutionary tracks. No strong constraint 
currently exists on the thick disc age. We assume an age 
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of 14 Gyr. An age of 11 Gyr, which is slightly older than 
the disc and younger than the halo, does not give signifi- 
cantly different results. 

The thick disc metallicity can be chosen between —0.4 
and —1.5 dex in the simulations. The standard value 
—0.7 dex is usually adopted, following in situ spectroscopic 
determination from Gilmore et aL (1995) and photometric 



We then compared the number of stars produced by 
the model with the observations in the selected region of 
the plane (magnitude, colour) and computed the likeli- 
hood that the observed data fits the model (following the 
method described in Bienayme et al., 1987a, appendix C). 



star count determinations (Robin et al., 1996; Buser et al 
1999|)~|The low metallicity tail of the thick disc seems to 



represent a weak contribution to general star counts (Mor- 
rison. Hgf ). it was neglected here. An internal metallic- 



ity dispersion among the thick disc population is allowed. 
The standard value for this dispersion is 0.25 dex with no 
metallicity gradient. 

The thick disc density law is assumed to be a truncated 
exponential: at large distances the law is exponential, at 
short distances it is a parabola ( Robin et al., 1996 ). This 
formula, given in equation 111 ensures the continuity and 
derivability of the density law (contrary to a true expo- 
nential) and eases the computation of the potential. 
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Three parameters define the density along the z axis : 
h z , the scale height, po the local density and xi the dis- 
tance above the plane where the density law becomes ex- 
ponential. This third parameter is fixed by continuity of 
p(z) and its derivative. It varies with the choice of scale 
height and local density following the potential. 



3.3. The spheroid 

We assume a homogeneous population of spheroid stars 
with a short period of star formation. We thus use the 
Bcrgbush fc VandcnBerg (199$ oxygen enhanced models, 
assuming an age of 14 Gyr (until more constraints on the 
age are available), a mean metallicity of —1.7 dex and 
a dispersion of 0.25 about this value. No galactocentric 
gradient is assumed. The colours are obtained from model 



atmospheres of Lejeune et al. (1997, 1998) 



The density law and IMF slope is the one deter- 
mined from deep star counts in numerous directions, as 



described in Robin et al. (2000). This is a power law 
dn/dm ex m -19 , an axis ratio of 0.7 and a local density 
of 1.64 10~ 4 stars pc~ 3 , excluding the white dwarfs. 



4. Data analysis method 

Population synthesis simulations have been computed in 
each observed field using photometric errors as close as 
possible to the true observational errors, generally grow- 
ing as a function of the magnitude and assumed to be 
Gaussian. Monte Carlo simulations were done in a solid 
angle larger or equal to the data in order to minimize the 
Poisson noise. 



The likelihood has been computed for a set of models, 
with varying thick disc parameters : scale height range 
400 to 1400 pc, scale length range 2 to 4 kpc and the IMF 
slope a from -0.25 to 2. In place of the local density we 
used a new parameter to try to overcome the degeneracy 
between the scale height and the local density. Since the 
number of stars is expected to vary as the volume times the 
scale height times the local density, we used the following 
density parameter, df : 
df = local density x scale height 2 

While this parameter does not reduce the degeneracy in 
shallow star counts (magnitude smaller than 20), it breaks 
the denegeracy in deeper counts, as can be seen in the 
likelihood contours in figure |^. The df parameter varies 
from to 2. 

The confidence limits of the estimated parameters are 
determined by the likelihood level which can be reached 
by random changing of the sample: a series of simulated 
random samples are produced using the set of model pa- 
rameters. The rms dispersion of the likelihood about the 
mean of this series gives an estimate of the likelihood fluc- 
tuations due to the random noise. It is then used to com- 
pute the confidence limit. Resulting errors are not strictly 
speaking standard errors; they give only an order of mag- 
nitude. 



5. Results 

5.1. Luminosity and mass function 

The IMF slope is best constrained when separately study- 
ing the fields, as the data do not cover the same mass range 
of thick disc stars. Shallow and deep fields give constraints 
on different parts of the luminosity function. Thick disc 
stars in DENIS fields have masses greater than 0.6M Q . 
The mass range of thick disc stars in deep counts is 0.2 
to 0.6M Q . The deepest field towards the North Galactic 
Pole (SA57) is dominated by stars with masses between 
0.2 and 0.4M Q . Figure || shows iso-contour likelihoods as a 
function of scale height h and density df for different IMF 
slopes, for DENIS fields, deep fields, and SA57, separately. 
An IMF slope a > 1.25 does not allow an acceptable solu- 
tion for all the fields. However, a lower IMF slope, a = 0.5, 
gives an agreement for all three magnitude intervals. 

Unlike DENIS fields, the best solution for deep fields 
is very sensitive to the IMF slope. This is also shown in 
figure U (solid lines). In a recent paper ( |Robin et al., 2000 ), 
we found an IMF slope a = 1.7 for the thick disc from 
analysis of a set of deep fields only. It appears now that 
a degeneracy between the parameters was not overcome 
and that the choice of a slightly larger scale height led to 
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a higher IMF slope (dotted lines in figure ||). However, we 
lacked a wide enough mass range to put a strong constraint 
on it. The solution given in the previous paper is not viable 
when we compare it to the DENIS star counts (dotted 
lines in figure ||). The only acceptable solution for the 
combination of deep and shallow star counts is a smaller 
IMF slope. 

We also tried to combine several IMF slopes in differ- 
ent mass intervals. The likelihood value is slightly better 
when a = 2 for DENIS fields and a = 1 for deep fields 
(see table 2). We tested an IMF with a change of slope 
at 0.6M©(My = 8), m~ 2 at higher masses and m _1 at 
lower masses, but the global likelihood value is smaller 
than when considering a single power law m" 5 (table 2) 
and the best solutions for the DENIS, deep, and Sa57 
fields are not in agreement. 

Paresce & Dc Marchi (2000) emphasize that the IMF 
of globular clusters of similar abundances, such as 47 Tuc, 
have the shape of a lognormai distribution: it rises as 
m~ 16 in the range 0.3-0.8M Q , then drops as m -0 2 be- 
low 0.3Mq. We find that such a lognormai distribution 
does not give a single solution acceptable for the DENIS, 
deep, and Sa57 fields (see also table 2 for the likelihood 
values) . 

The fainter thick disc stars in the deep fields are in 
the mass range 0.2-0.5M Q . The deepest bin 22-24 in Sa57 
contains thick disc stars with masses from 0.1 to 0.4M Q . 
This field alone does not give enough constraints to deter- 
mine if an IMF with a change of slope around 0.3M Q would 
give a better agreement because of the sample size. Large 
scale surveys at this depth, like MEGACAM or VISTA 
projects, would be necessary to definitely choose between 
several power laws or a lognormai IMF. 



5.2. Scale height and scale length 

As shown in figure |], our best fit model for all the fields 
together gives a scale height of 800 pc with df = 1. Iso- 
contours at 1, 2 and 3 a are also plotted. Assuming this 
scale height and density, the maximum likelihood is ob- 
tained for a scale length of 2500 pc, but smaller or greater 
values are still acceptable. We now consider separately the 
fields in the anticentre (135° < 1 < 225°) and the center (1 
> 315° or 1< 45°), at medium galactic latitude (|b| < 50°) 
where the effects of the scale length are most important. 
Iso-contours at 1 a are plotted in figure || with solid lines 
for the centre fields and dashed lines for the anticentre 
fields, the scale length ranging from 2000 pc to 3500 pc. A 
scale length of 2500 pc gives the best agreement between 
the centre and anticentre fields. Disc scale length is about 



5.3. Sensitivity to other parameters 

5.3.1. The Age of the thick disc 

Considering a younger thick disc of 11 Gyr instead of 14 
Gyr gives the same results over the scale height, density 
and scale length. The IMF slope that allows us to recon- 
cile deep and shallow fields is slightly different, a = 0.75, 
but still within our error bars. However, the absolute like- 
lihoods are not as good as the ones obtained with an age 
of 14 Gyr. This parameter should mainly be determined 
from the turnoff position with accurate enough counts in 
a homogeneous system. The color shift in I-J for a 8 Gyr 
to a 14 Gyr thick disc is only of 0.06 magnitude at the 
turnoff, out of reach at the DENIS precision. 

5.3.2. The spheroid IMF 

We have consid ered an IMF slope a — 1.9 for the spheroid, 
as derived by Robin et al. (2000 ) from the deep star 
counts. Even con sidering an IMF slo pe as low as a — 0.75 
for the spheroid ( Gould et al., 1998 ) does not change our 
results concerning the thick disc IMF slope. The best fit 
parameters are h = 850 pc and df = 1.1, well within our 
error bars. 



5.3.3. Binary effects 

The thin disc luminosity function that we used is a nearby 
luminosity function that does not take into account the 
fact that we may observe systems instead of single stars. 
Kroupa (2000b] ) showed that the difference between the 
nearby luminosity function and the system one comes from 
the binary fraction, nearby systems being computed as 
true single stars, while systems are not resolved in remote 
star counts. Hence, a correction is to be applied to the 
nearby luminosity function to take into account this effect 
in the counts. Following Kroupa, we applied a correction 
to the luminosity function used by the model in order to 
match the system luminosity function, as shown in fig- 
ure ?? (upper curves). Only the deepest fields could be 
sensitive to this effect. When we apply this correction, the 
best fit thick disc parameters remain unchanged. 

For the thick disc, the luminosity function considered 
in the model is the luminosity function of single stars. 
Knowing little about the binarity in the thick disc, we tem- 
porarily applied the same correction as for the thin disc 
on the luminosity function (see figure ??, lower curves). 
The best fit is obtained for the same density parameters 
h and df, but it slightly displaces the best IMF slope to 
a = 0.75. The likelihood values are slightly higher than 
those obtained with a single luminosity function (table 2) . 



2500 pc, from the most recent studies (Robin et al., 1992; 



Fux fc Martinet, 1994; Ruphy et al., 1996). The thick disc 



one seems to be very similar. The lack of deeper fields close 
to the anticentre does not allow us to better constrain this 
parameter or to reveal the presence of a flare. 



6. Discussions and conclusions 

We have estimated the thick disc density law parameters 
and mass function by using a wide set of data at high and 
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intermediate galactic latitudes, in the visible and the near 
infrared. The best fit model has a scale height of 800 pc, a 
scale length of 2500 pc and a density of 10 -3 stars pc -3 or 



galactic field IMF for the disc has a slope of 1.3±0.5 at 
O.lM Q <m<O.5M and 2.3±0.3 at O.5M <m<lM^ 



7.1 IP" 4 M pc~ 3 for Mv < 8, that is 6.2% of the thin disc 



density. This result confirms the values obtained in 1996 
from a smaller number of fields and mass ranges ( |Robirj 
et al, jUjj ). 

For the first time, we determined the thick disc mass 
function over a large mass range by the study of shallow 
star counts dominated by stars at the turnoff combined 
with deep star counts. The IMF of the thick disc seems 
to follow a power law dn/dm oc m~ in the mass range 
0.2-0.8M Q . We found no evidence of a change of slope at 
lower masses, but we only have one field deep enough to 
constrain the IMF at low masses. 

The only point of comparison for the thick disc mass 
function is the globular clusters of similar metallicities, 
where the conditions of star formation could be compara- 
ble. However, clusters are subject to bias with regard to 
the field because of possible mass segregation effects. From 
seven clusters, Piotto & Zoccali (1999) found a mean a of 



0.89 at masses below 0.6M Q but the range covers 1.22 to 
0.53, which is compatible with our measurements in the 
field. 

Paresce & De Marchi (2000) found for similar clusters 
an indication of a lognormal IMF rather than a power law. 
If fitted by a power law in the mass range 0.3-0.8M©, their 
IMF should be approximated by cv=1.6, but going down 
to a=0.2 at m<0.3 M Q . The result clearly depends on the 
way an IMF, which globally is not a power law, is fitted in 
different mass intervals by portions of power laws. While 
our result seems not to favor a lognormal IMF compared 
to a power law slope in the mass range 0.2-0.8M Q , more 
data at lower masses may change our conclusion in the 
future. 

Values of the spheroid IMF slope in the field range 



between 1.7±0.3 from a small local sample (Chabrier & 



Mera, 1997j ), 0.75 from HST star counts ( pould ct al 
1998| ), and the higher value a=1.9 determined from wide 



field star counts (Robin ct al., 2000). However, the latter 
is valid at masses m>O.3M while the former go deeper 
to masses of about 0.1M Q and are less sensitive to higher 
masses because of the small size of the samples. A change 
of slope between intermediate masses and low masses, as 
found in the disc, may explain the apparent discrepancy, 
which is also partly due to the Poisson noise in these rel- 
atively small samples. 

Scalo (1998) suggests it is not valid to rely upon a 
mean IMF, since variations are too high from one mea- 
surement to the other at the present time, but he iden- 
tifies no tendency correlated with any physical parame- 
ter. On the other hand, Kroupa (2000a| ) finds that young 
clusters seem to have a steeper slope for m< lM Q and an- 
cient globular clusters have a > but closer to than 
the field. He concludes that star formation produces rel- 
atively more low mass stars at later galactic epochs. His 



If we consider the alpha-plot (figure 14 of Kroupa 
(2000a)), our measurement of a thick disc mass function 



similar to globular clusters seems to corroborate Kroupa's 
conclusion. The thick disc mass function at low masses 
seems to be flatter than the thin disc and comparable with 
the spheroid. Figure fj] shows the alpha-plot versus metal- 
licity for low-mass stars (m<3M Q ), for which the effect 
is clearer. However, when taking into account other halo 



IMFs measured in the field, like Chabrier & Mera (1997), 



the correlation is less clear because of larger error bars due 
to the small size of the sample, as seen in the alpha-plot 
versus mass (figure [7]). Clearly, the answer depends on a 
better determination of the IMF in the field, which can be 
done by exploring a larger mass range (going deeper) and 
having larger samples. This is promising in that several 
large scale surveys are on the way or planned for the near 
future. 
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Fig. 3. Colour distributions of a DMS field in the magni- 
tude bin 1=20-22 and a DENIS field in the magnitude bin 
1=14-16. The dots show the observations. Thick lines are 
predicted number of stars by the model assuming a = 0.5, 
h = 800 pc and df = 1. Thin lines: a = 2, h = 800 pc, df 
= 1. Dotted lines: a = 1.7, h = 1100 pc, df = 0.7. 
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Fig. 4. Iso-contour likelihoods at 1, 2, and 3 ct as a func- 
tion of scale height h and density df for all the fields to- 
gether, considering an IMF slope a = 0.5. 
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Table 1. Characteristics of the selected fields. 



Reference 


Field 




Area 


Bands 


Ma 


gnitude 




coordinates 


(deg 2 ) 




ran 


ge 


Ojha et al. (1994a,b, 1996) 


1=4°, b=+47° 


15.5 


B,V 


V= 


15-17 




1=278°, b= 


=+47° 


20.8 


B,V 


V= 


16-18.5 




1=168°, b= 


=+48° 


7.13 


B,V 


V= 


15-18 


Chiu (1980) 


1=189°, b= 


=+21° 


0.10 


B,V 


v= 


18-20 


Yee et al. (2000) 


1=166°, b= 


=-55° 


0.39 


v,i 


v= 


17-20 


Borra & Lepage (1986) 


1=197°, b= 


=+38° 


0.25 


B,V 


v= 


20-22 


Reid & Majewski (1993) 


North Galactic Pole 0.30 


B,V 


v= 


19-22 


DMS 


1=129°, b= 


=-63° 


0.14 


v,i 


v= 


18-22 




1=248°, b= 


=+47° 


0.08 


v,i 


v= 


18-22 




1=337°, b= 


=+57° 


0.16 


v,i 


v= 


18-22 




1=77°, b= 


+35° 


0.15 


v,i 


v= 


18-22 




1=52°, b= 


-39° 


0.15 


v,i 


v= 


18-22 




1=68°, b= 


-51° 


0.15 


v,i 


v= 


18-22 


Bouvier et al. (1998) 


1=117°, b= 


=-23° 


1.48 


I,R 


R= 


17-22.4 




1=206°, b= 


=+32° 


0.55 


I,R 


R= 


17-22 


Sa57 


1=69°, b= 


+85° 


0.16 


v,i 


v= 


20-24 



Table 2. Likelihood values obtained for the best fit parameters h and df when considering different thick disc IMF. 
The best fit parameters h and df are shown in figure | for the DENIS, deep and Sa57 fields in the case of an IMF 
following a single power law. 



dn/dm 


DENIS fields 


deep field 


SA57 


all fields 


m-°- 5 


-1830+200 


-2500+150 


-92+10 


-5300+400 


m-°' 75 


-1790 


-2420 


-88 


-5250 


m" 1 


-1760 


-2370 


-85 


-5370 


m" 1 - 5 


-1730 


-2370 


-79 


-6260 


m" 2 


-1730 


-2480 


-77 


-8130 


m" 2 (m>O.6M ), m" 1 (m<O.6M ) 


-1730 


-2570 


-85 


-7230 


i [log(m/o.34)] 2 , 
m ex P 2(0.4)* (lognormal) 


-1710 


-2510 


-89 


-6310 


m -o.75 _|_ ]3j naX y correction 


-1820 


-2340 


-87 


-5100 
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Galactic longitude 1 



Fig. 1. Galactic coordinates of the selected fields. Characteristics of the fields plotted with dots are given in table 1. 
The lines represent parts of DENIS strips. 
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Fig. 5. Iso-contour likelihoods at 1 a as a function of scale height h and density df for different scale lengths 1. Plus 
and solid line: fields towards the galactic centre. Diamond and dashed line: fields in the anticentre direction. 
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Fig. 7. IMF slope a versus mass in MQand metallicity. 7, T and x are also power-law indices that may be used by 



other authors. Thin lines and filled symbols: globular clusters (triangles down: Paresce & De Marchi (2000), triangles 
up: Piotto & Zoccali (1999)). Thick lines and filled symbols: spheroid field stars (circles: Robin et al. (2000^7 diamonds : 
Chabrier & Mera (1997), squares: Gould et al. (1998| )). plus: thick disc stars. Open circles: disc stars flKroupa, 2001 ). 



